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ABSTRACT 



We present the extension of our NextGen model atmosphere grid to the 
regime of giant stars. The input physics of the models presented here is nearly 
identical to the NextGen dwarf atmosphere models, however spherical geometry 
is used self-consistently in the model calculations (including the radiative 
transfer). We re- visit the discussion of the effects of spherical geometry on the 
structure of the atmospheres and the emitted spectra and discuss the results of 
NLTE calculations for a few selected models. 
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1. Introduction 

In a recent paper ( |Hauschildt, Allard, fc Baron, 199"9| , hereafter: NG) we have presented 
a grid of plane parallel model atmospheres for dwarf stars (log(g) > 3.5) and in the 
temperature range 3000 < T c r < 10000 K. The NG models are intended for the analysis 
and modeling of main sequence stars and, where applicable, sub-giants. Spherical extension 
effects become more important for smaller gravities, which is the reason why the NG-grid 
stopped at log(g) = 3.5. In this paper, we present an extension of the NG grid to smaller 
gravities by including the effects of spherical geometry in the model calculations. This 
includes, for example, the calculation of the hydrostatic structure in spherical geometry and 
the effects of spherical radiative transfer on the model atmospheres. 

In order to retain compatibility with the NG grid, the NG-giant grid presented here 
uses the same parameterization of the element abundances as the NG grid. Therefore, the 
models presented here are not applicable to chemically peculiar stars such as, e.g., Carbon 
stars. However, there are many possible applications for these models. For example, 
the Period-luminosity (PL) relation of Cepheids is of fundamental importance for the 
determination of the extragalactic distance scale. [Alibert et al. (1999| ) have used the 



atmosphere models presented in this paper to connect evolutionary models to observed 
colors for Cepheids with very good results. 

In the next section we give a brief overview over the model construction and the 
differences from the NG grid. Then we discuss the results, in particular the effects of 
spherical symmetry and we end with a summary of the paper. 



2. Model calculations 



We have calculated the models presented in this paper using our multipurpose model 
atmosphere code PHOENIX, version 10.5. Details of the code and the general input physics 
setup are discussed in |Hauschildt, Allard, fc Baron ( 1999]) and |Hauschildt fc Baron (1999|) 
and references therein. The model atmospheres for Cepheids presented here were calculated 
with the same general input physics. However, they use spherical geometry (including 
spherically symmetric radiative transfer) rather than plane parallel geometry. For giant 
models with low gravities (log(g) < 3.5), this can be important for the correct calculation 



of the structure of the model atmosphere and the synthetic spectrum (c.f., e.g. |Aufdenberg 



et al., 1998a, b, for a discussion of spherical effects in B giants) 



Our combined molecular line list includes about 500 million molecular lines. These 
lines are treated with a direct opacity sampling technique where each line has its individual 
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Voigt (for strong lines) or Gauss (weak lines) line profile. The lines are selected for every 
model from the master line list at the beginning of each model iteration to account for 
changes in the model structure (see [Hauschildt, Allard, fc Baron, 199§|, or details of the line 



selection process). This procedure selects about 190 million molecular lines for a typical 
giant model atmosphere with T e g « 3000 K. Accordingly, we generally use the parallelized 
version of PHOENIX ( |Hauschildt, Baron, fc Allard, 19971; Baron fc Hauschildt, 1998|) to 



perform calculations efficiently on parallel supercomputers. Details of the calculations are 
given in the above references and are not repeated here. 

The parameterization of giant models requires an additional parameter compared to 
plane parallel model atmospheres. This complicates the task of relating theoretical models 
to observed data, see [Baschek, Scholz, fc Wehrse (1991| ) for a discussion of these issues. In 
the models presented here, we set the stellar radius R by the condition g gra , v = 

GM/R 2 , 

where we define g grav as the gravitational acceleration at r stc j = 1, G is the constant of 
gravity, and r st( j is the optical depth in the continuum at 1.2 fi. The luminosity L of the 
model is then given by L = ATiR 2 aT e ^. For convenience, our model grid is based on the 
set of parameters ( T cS , log(g), M, [M/H]). The above formulae and the structures of the 
model atmospheres can be used to transform them to any target set of parameters (e.g., for 
a different definition of r st d)- 



2.1. Temperature correction procedure 

We iterate for the temperature structure of the atmosphere using a generalization of 
the Unsold-Lucy temperature correction scheme to spherical geometry and NLTE model 
calculations. This has proven to work very well even in extreme NLTE cases such as nova 
and supernova atmospheres. The temperature correction procedure also requires virtually 
no memory and CPU time overhead. The Unsold-Lucy correction scheme (see |Mihalas,' 



1970| , for a discussion of this and other temperature correction schemes), uses the global 
constraint equation of energy conservation to find corrections to the temperature that 
will fulfill energy conservation better than the previous estimate. We have found it to be 
more stable than a Newton-Raphson linearization scheme and it allows us to separate the 
temperature corrections from the statistical equations. The latter property is extremely 
useful in NLTE calculations where it allows us to use a nested iteration scheme for the 
energy conservation and the statistical equilibrium equations. This significantly increases 
numerical stability and allows us to use much larger model atoms than more conventional 
methods. 



To derive the Unsold-Lucy correction, one uses the fact that the ratios of the wavelength 
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averaged absorption and extinction coefficients 

k p = n x Bid\\ /B 

Xj 



(J^KxJidXj/J (2) 
J^XxFtdX^/F (3) 



(4) 

(where B, J,F denote the wavelength integrated Planck function, mean intensity and 
radiation flux, respectively) depend much less on values of the independent variables than 
do the averages themselves. 

Dropping terms of order (v/c), one can then use the angular moments of the radiative 
transfer equation to show that in order to obtain radiation equilibrium B should be 
corrected by an amount 

SB(r) = —{kjJ - k p B + S/(4tt)} (5) 

Kp 

-HL{2(H(t = 0)-H (t = 0)) 

Kp 



~~n*J Rq r ' 2 XF {H{r ' ] ~ Ho{r,)) dr ' h 




fq 

where H = F/Atc, H (t) is the value of the target luminosity at that particular depth point 
(variable due to the velocity terms in comoving frame radiative transfer calculations and 
non-mechanical energy sources, the total observed luminosity Hq(0) is an input parameter), 
Here, q is the "sphericity factor" given by 

1 

q = exp 

where r core is the inner radius of the atmosphere, R is the total radius, /(r) = K(t)/J(t) 
is the "Eddington factor" , and K = J fi 2 I dp, is the second angular moment of the mean 
intensity. S describes all additional sources of energy such as mechanical energy supplied 
by winds or non-thermal ionization due to 7-ray deposition. 

The first term in Eq. || corresponds simply to a A iteration term and will thus provide 
temperature corrections that are smaller than required in the inner parts of the atmosphere, 
but will be accurate in the outer, optically thin parts. The second term of Eq. [5], however, 
is the dominant term in the inner parts of the atmosphere. It provides a very good 
approximation to the temperature corrections AT deep inside the atmosphere. Following 
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Unsold (1968 ), we have found that it is sometimes better to modify this general scheme by 
excluding the contributions of extremely strong lines in the opacity averages used for the 
calculations of AT because they tend to dominate the average opacity but do not contribute 
as much to the total error in the energy conservation constraint. 



3. Results 

We have calculated a grid of solar abundance (table 5 of Jaschck fc Jaschck, 1995) 



model atmospheres. In addition, we have calculated grids with metallicities [M/H] = —0.3, 
—0.5 and —0.7 to facilitate comparison to observed LMC and SMC giants. The models span 
a range of 2000 K < T e g < 6800 K and 0.0 < log(g) < 3.5. The models are not available for 
all combinations of these parameters, e.g., at high T e s and low log(g) the total radiative 
acceleration is larger than the local gravity and hydrostatic models become invalid. In order 
to minimize CPU time, the models were calculated for a fixed mass of M — 5M . In the 
following, we discuss the effects of different model parameters on the resulting structures 
and spectra. The effects of the model parameters on the colors are discussed in the context 
of evolutionary models by |Alibert et al. (1999| ). 



3.1. Structure of the atmospheres 

If Fig. [I] we show the temperature structure for models with varying log(g) for 3 
effective temperatures. For the models with T e g = 5600 K, and the smallest gravity, the 
atmosphere does not become convective within the maximum optical depth of the model, 
r std — 100. The apparent flattening of the temperature drop in the outer layers of the 
atmosphere is a result of the automatic "compression" of the optical depth scale toward the 
outer edge of the model atmosphere. This means that the radii of the different optical depth 
points are much closer to each other in the outer atmosphere compared to r s td ~ 1- This 
is most pronounced in the coolest models shown, whereas the "scale compression" is nearly 
non-existent for the model with T e g = 5600 K. We have also compared model atmospheres 
calculated with the same parameters (T e g, log(g)) but for different stellar masses. For 
T e s = 3600 K, log(g) = 0.0 and solar abundances, we find only small temperature differences 
for stellar masses between 2.5 and 7.5 M Q . However, the absolute differences in the synthetic 
spectra are equivalent to a significant change in log(g) (see below) whereas the relative 
differences in the spectra are small. The radial extension of the atmospheres, here defined 
simply as the ratio of the outer radius R out and the inner radius Ri n of the atmosphere, is 
typically less than 20% (see Fig. §). We have ignored all models that become unstable due 
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to large radiative accelerations (low \og(g) models with very low or high temperatures), such 
atmospheres have to be described by stellar wind models (Aufdenberg et al., in preparation) 
and typically have much larger radial extensions. 



3.2. Sphericity effects 

It is important to assess the relevance of the effects of spherical geometry and radiative 
transfer on the synthetic spectra. For hot stars, Aufdenberg et al. (1998a ^rj|) have shown 



that a combination of geometric effects and line blanketing can resolve a long standing 
problem of discrepancies between synthetic and observed spectra. The effects of spherical 
symmetry on cool giant atmospheres has been investigated previously by [Scholz fc Tsuji| 
I984D for giant M and C stars, by piez, Brett, fc Nordlund (1992|) for M giants, and by 



J0rgensen, Johnson, fc Nordlund (1992|) for Carbon stars. Our model grid extends this 



work toward warmer stars and reduces the gap to the stars considered in |Aufdenberg et aT 



1998a|) . In addition, improved molecular line lists are now available and since |Aufdenbcrg| 



et al. (1998a| ) found a close connection between line blanketing and sphericity for hot stars, 



the improved molecular line data could affect the results. 

In Fig. ^ we compare normalized (to the same area) low-resolution (20A) spectra for 
3 sets of models. Each panel displays a NG-giant synthetic spectrum (full lines) and a 
synthetic spectrum calculated with the same parameters and input physics, but assuming 
plane parallel geometry. In this figure, only the model with T e s = 4000 K shows visible 
differences between the plane parallel and spherical cases. However, plotting the relative 
differences of the spectra in Fig. ^ reveals systematic differences even for low resolution 
spectra. For T e fj = 3000 K (top panel), we find that the spherical model emits more flux 
than the plane parallel model in the short wavelength regime. There are also some changes 
in the slopes of the spectra between 0.45 and 0.55 /im, which are hardly noticeable in the 
low resolution spectra. The middle panel shows the results for T e g = 4000 K. This model 
displays a different systematic change in the blue spectral region, the long wavelength part 
of the displayed wavelength interval does not reveal large systematic effects. The highest 
effective temperature shown, T e s = 5600 K, shows the smallest systematic changes. 

Figure |5| is similar to the previous plot but it shows the changes for a resolution 
of 2A (10 times better). From this figure it becomes clear that the changes are due to 
the effect of line blanketing: line overlap and individual line strengths change between 
spherical and plane parallel models. In low resolution spectra this produces a change in the 
pseudo-continuum that is formed by the overlapping spectral lines. In the T e g = 3000 K 
model the TiO lines are significantly affected. This shows that spectral analyses based on 
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high-resolution spectra need to account for spherical effects. 

For longer wavelengths we find that the differences between spherical and plane parallel 
models are smaller. Figures ^ and [7| demonstrate this for low resolution and medium 
resolution spectra, respectively. The changes are mostly due to the sensitivity of the CO 
lines to changes in the structure of the atmosphere, therefore, they are larger at lower 
effective temperatures. The figure with the higher resolution shows that there are significant 
changes in individual water lines for low T e g. At higher effective temperatures, individual 
atomic lines are up to 15% stronger in the plane parallel models (less emitted flux, since 
these are absorption lines). The conclusion is that high- resolution spectral analysis requires 
spherical models, whereas low resolution spectra or colors can be interpreted using simple 
plane parallel models. 

In order to demonstrate the reason for these differences, we display in Fig. § the 
electron temperatures as functions of our standard optical depth r st d, which is simply the 
optical depth in the continuum (b-f and f-f processes) at a reference wavelength of 1.2 /xm. 
The temperature structures for the 2 models with T e g- = 5600 K are nearly identical, the 
differences are less than about 80 K everywhere. In the model with the lowest effective 
temperature (3000 K) the differences are larger, close to 200 K at r st d ~ 10 -4 . This 
amounts to nearly 10% change in the absolute temperature, which is quite significant for 
the formation of the spectrum. Note that the line optical depth in the blue spectral region 
can be orders of magnitude larger than r st d and many "blue lines" actually form around 
Tstd ~ 10~ 4 . The model with T e g = 4000 K lies in between the two extremes and is not 
shown in the figure. The temperature structure is only a part of the full structure of the 
atmosphere, gas pressures and partial pressures change accordingly. For larger gravities, 
the differences between plane parallel and spherical models become smaller and for gravities 
log(g) > 3.0 the differences are small enough to be negligible for most applications. 

3.3. Spectra 

The changes of the synthetic spectra with gravity are illustrated in Fig. [| The 
differences are larger for the cooler model in the sense that at lower T e Q, higher gravity 
model atmospheres emit more optical flux than lower gravity models. This is caused by 
the enhanced near-IR absorption of water vapor in the models with larger gravity and 
hence larger gas pressures in the line forming region. For larger effective temperatures 
the differences diminish and are confined to individual gravity sensitive features in the 
spectrum. 
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Changes in the metallicity [M/H] result in changes of the synthetic spectra that depend 



strongly on the effective temperature, as shown in Fig. [It]. For high effective temperatures, 
the changes are small in the [M/H] range we have investigated. At lower T e g, however, the 
changes are dramatic as illustrated in the top panel of Fig. [TO. This is due to the increasing 



importance of molecules at lower temperatures. The concentration of molecules and thus 
their opacities depend strongly on the metal abundances. Molecules are less important at 
higher T e g and thus the spectra are less sensitive to metallicity changes. 

The very small sensitivity of the atmospheric structure on the mass of the star is 
mirrored by only small changes in the resulting low- resolution spectra (Fig. [11]). The 



differences correspond to scaling factors close to unity over a large wavelength range and 
are thus basically negligible for most purposes. The absolute changes in limited wavelength 
intervals of the spectra, however, can be equivalent to changes in log(g), thus the mass of 
the star needs to be considered as a parameters for applications that use absolute spectra. 



3.4. NLTE effects 

We have calculated a small number of NLTE models in order to investigate the 
importance of NLTE effects on the structure of the model atmospheres. The results for 



cooler models were discussed in Hauschildt et al. (1997|) and are not repeated here. Figures 



12| and [13] show an overview of selected NLTE species for models with T e g = 4000 K and 
5600 K for log(g) = 0.0 and solar abundances. The total number of NLTE levels in each 
model is 4532 with a total of 47993 primary NLTE lines (see |Hauschildt fe Baron, 1999|, and 



references therein for details). The following species (and number of levels) were treated in 
NLTE: H 1 (30), Mg I (273), Mg 11 (72), Ca 1 (194), Ca 11 (87), Fe I (494), Fe II (617), 
O I (36), O II (171), Ti I (395), Ti II (204), C I (228), C II (85), N I (252), N II (152), 
Si I (329), Si II (93), S I (146), S II (84), Al I (111), Al II (188), K I (73), K II (22), Na I 
(53), and Na II (35). For most of the species, the departure coefficients are always close to 
unity, in particular for species with resonance lines and photoionization edges in the UV 
part of the spectrum. The species shown in Figure [12] and [13] are the ones with the most 
pronounced departures from LTE. The departures are generally too small to significantly 
affect the structure of the atmospheres. However, NLTE effects do change the profiles of 
individual lines as shown in |Hauschildt et al. (1997| ) for Ti I lines. Therefore, abundance 



analyses of individual elements should take NLTE effects into account whenever possible. 
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4. Summary and Conclusions 

In this paper we presented the extension of our grid of dwarf model atmospheres 
( |Hauschildt, Allard, fc Baron, 199~P| ) to giant stars. The models were computed using 
spherical symmetry, including spherical radiative transfer. We discuss the results of 
the calculations. We have calculated a small number of NLTE models and verified the 



conclusions of [Hauschildt et al. (1997|) : In the considered range of effective temperatures, 
NLTE effects do not alter the structure of the atmospheres significantly but they can 
change the profiles of individual lines. NLTE should therefore be considered in abundances 



analyses. Our grid has been used successfully in evolutionary models of Cepheids flAlibert 



|et al., 1999 ). We provide the model structures and spectra through the WWW and 



anonymous FTP for general use, see [http : / /dilbert . physast . uga . edu/~| yeti or 



ftp : //calvin . physast . uga . edu/pub/NG-giant|. In the next version of the grid we will 



use improved molecular line opacities as well as dust opacities for low effective temperatures. 
In addition, stellar winds models are needed for those models that are unstable against 
radiation pressure, which occurs at both high and low effective temperatures. 
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5. Figures 



Fig. 1. — Temperature structures for solar abundance models with T cS = 5600 K, 4200 K, 
and 3000 K for varying gravities (as indicated). r std is the optical depth in the continuum 
(b-f and f-f processes) at a wavelength of 1.2 /im. 
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Fig. 2. — Radial extensions R ut/Rm for solar abundance models with T c g = 5600 K, 4200 K, 
and 3000 K (as indicated). 

Fig. 3.— Comparison between solar abundance giant models calculated using spherical 
geometry and radiative transfer (full curves) and plane parallel geometry and radiative 
transfer (dotted curves) in the blue/optical spectral region. The resolution of the spectra 
has been reduced to 20A. 

Fig. 4. — Relative flux change between spherical and plane parallel model calculations. The 
y-axis shows (f p — f s )/f s , where f p is the flux of the plane parallel model and f s is the flux 
calculated for the spherical model. The resolution of the spectra was reduced to 20A. 

Fig. 5. — Relative flux change between spherical and plane parallel model calculations. The 
y-axis shows (f p — / s )// s , where f p is the flux of the plane parallel model and f s is the flux 
calculated for the spherical model. The resolution of the spectra is 2A. 

Fig. 6. — Relative flux change between spherical and plane parallel model calculations. The 
y-axis shows (f p — f s )/fs, where f p is the flux of the plane parallel model and f s is the flux 
calculated for the spherical model. The resolution of the spectra was reduced to 20A. 

Fig. 7. — Relative flux change between spherical and plane parallel model calculations. The 
y-axis shows (f p — f s )/ f s , where f p is the flux of the plane parallel model and f s is the flux 
calculated for the spherical model. The resolution of the spectra is 2A. 

Fig. 8. — Temperature structures for spherical (full curves) and plane parallel (dotted lines) 
model calculations. r std is the optical depth in the continuum (b-f and f-f processes) at a 
wavelength of 1.2 /im. 

Fig. 9. — Sensitivity of the synthetic spectra to gravity changes for solar abundance models 
with T cS = 5600 K, 4200 K, and 3000 K (as indicated). The resolution of the spectra has 
been reduced to 20A. 

Fig. 10.— Sensitivity of the synthetic spectra to metallicity changes for models with 
T cS = 5600 K, 4200 K, and 3200 K (as indicated) and \og(g) = 0.0. For the T eff = 3200 K 
model the metallicities [M/H] = 0.0 and [M/H] = —0.3 are shown whereas for the hotter 
two models the metallicities shown are [M/H] = 0.0 and [M/H] = —0.7. The resolution of 
the spectra has been reduced to 20A. 



Fig. 11. — Synthetic spectra at 20A resolution for models with T eS = 3600 K, log(g) = 0.0 
solar abundances, and stellar masses of 7.5 M Q and 2.5 M . 



Fig. 12. — Overview over selected departure coefficients for a NLTE model with T eS 
4000 K, log(g) = 0.0, and solar abundances. 



Fig. 13. — Overview over selected departure coefficients for a NLTE model with T e g 
5600 K, log(g) = 0.0, and solar abundances. 



